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Abstract 

Periodic nonuniform sampling is a known method to sample spectrally sparse signals below the 
Nyquist rate. This strategy relies on the implicit assumption that the individual samplers are exposed 
to the entire frequency range. This assumption becomes impractical for wideband sparse signals. The 
current paper proposes an alternative sampling stage that does not require a full-band front end. Instead, 
signals are captured with an analog front end that consists of a bank of multipliers and lowpass filters 
whose cutoff is much lower than the Nyquist rate. The problem of recovering the original signal from 
the low-rate samples can be studied within the framework of compressive sampling. An appropriate 
parameter selection ensures that the samples uniquely determine the analog input. Moreover, the analog 
input can be stably reconstructed with digital algorithms. Numerical experiments support the theoretical 
analysis. 

I. Introduction 

Radio frequency (RF) technology enables the modulation of a narrowband signal by a high carrier 
frequency. As a consequence, manmade radio signals are often sparse. That is, they consist of relatively 
small number of narrowband transmissions spread across a wide territory of spectrum. A convenient 
description for these signals is the multiband model where the frequency support of a signal resides 
within several continuous intervals in a wide spectrum but vanishes elsewhere. 
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It has become prohibitive to sample modern multiband signals because their Nyquist rates may exceed 
the specifications of the best analog-to-digital converters (ADCs) by orders of magnitude. As a result, 
any attempt to acquire a multiband signal must exploit its structure in an intelligent way. 

Previous work on multiband signals has shown that it is possible to reduce the sampling rate by 
acquiring samples from a periodic but nonuniform grid [1]. Multi-coset sampUng, a specific strategy of 
this type, was analyzed in [2], which established that exact recovery is possible when the band locations 
are known. The bUnd case, in which the band locations are unknown, has been extensively studied in 
[3]. Unfortunately, the sampling front ends proposed in [l]-[3] are impractical for wideband applications 
because they require ADCs whose sampling rate is matched to the Nyquist rate of the input signal, even 
when the average sampling rate is much lower. Other limitations are described in Section II-B. Another 
recent work [4] has partially overcome these shortcomings using a hybrid optic-electronic system at the 
expense of size and cost. 

In this paper, we analyze a practical sampling system inspired by the recent work on the random 
demodulator [5]. This system multiplies the input signal by a random square wave alternating at the 
Nyquist rate, then it performs lowpass filtering, and samples the signal at a lower rate. Our system 
consists of a bank of random demodulators miming in parallel. We show that, for an appropriate choice 
of parameters, our system uniquely and stably determines a multiband input signal. Moreover, we describe 
digital algorithms for reconstructing the signals from the parallel samples. 

We continue with an outline of the paper. Section II reviews essential background material. In Sec- 
tion III, we describe the system design and a frequency-domain analysis that leads to an infinite measure- 
ment vectors (IMV) system. Applying ideas from [6], we reduce the problem of locating the frequency 
bands to a finite-dimensional compressive sampling problem. We then derive an appropriate choice of 
parameters for the sampling system. Section IV presents our numerical experiments, which demonstrate 
that the system permits stable signal recovery in the presence of noise. 

II. Formulation and Background 

A. Design Goals for Efficient Sampling 

Let x{t) be a real-valued, finite-energy, continuous-time signal, and let X{f) = x{t) exp(— j27r/t)di 
be its Fourier transform. We treat a multiband signal model M in which x{t) is bandlimited to = 
[— /nyq/2, /nyq/2] and the support of X{f) consists of at most 2N frequency interval whose widths do 
not exceed B. Fig. I depicts a typical communication appHcation that obeys this signal model. 



Fig. 1: Three RF transmissions witli different carriers fi. The receiver sees a multiband signal (bottom drawing). In this 
example N = 3, and the modulation techniques of the transmitters determine the maximal expected width B. 



We wish to design a sampling system for signals from our model M. that satisfies the following 
properties: 

1) The sampling rate should be as low as possible; 

2) the system has no prior knowledge of band locations; and 

3) the system can be implemented with existing devices. 
We will call this type of sampling stage efficient. 

The set is a union of subspaces corresponding to all possible signal supports. Every x{t) e M hes 
in one of these subspaces. Detecting the exact subspace, prior to samphng, may be impossible or too 
expensive to implement. An efficient system should therefore be blind, in the sense that band locations 
are not assumed to be known. 

The lowest (average) samphng rate that allows bhnd perfect reconstruction for all signals in is ANB 
samples/sec [3]. This rate is proportional to the effective bandwidth of x{t), and it is typically far less 
than the Nyquist rate /nyq> which depends only on the maximum frequency in x(t). See Section HI-C 
for more discussion. 

Our previous work describes blind reconstruction of x{t) £ M from multi-coset samples taken at the 
minimal rate [3]. The next section details the practical hmitations of the multi-coset strategy, which make 
it inefficient for wideband signals. 

B. Practical Limitations of Multi-Coset Sampling 

Multi-coset sampling involves periodic nonuniform samphng of the Nyquist-rate sequence x(nrNYQ), 
where Tnyq = 1//nyq- The ith coset takes the i\h value in every block of L consecutive samples. 
Retaining only p < L cosets, indexed by C = {ci}^^-^^, gives p sequences 




otherwise. 



n = mL + Ci,m G 



(1) 
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with an average sampling rate p/(LTnyq), which is lower than the Nyquist rate. 

To explain the practical limitations of this strategy, we observe that standard ADC devices have a 
specified maximal rate r, and manufactures require a preceding low-pass filter with cutoff r/2. Distortions 
occur if the anti-aliasing filter is not used, since the design is tailored to r/2-bandlimited signals and 
has an internal parasitic response to frequencies above r/2. To avoid these distortions, an ADC with r 
matching the Nyquist rate of the input signal must be used, even if the actual sampUng rate is below 
the maximal conversion rate r. In multi-coset sampUng, each sequence Xcjn] corresponds to uniform 
sampling at rate 1/(LTnyq), whereas the input x{t) contains frequencies up to /^fYQ/2. Acquiring Xcjn] 
is only possible using an ADC with r = /nyq, which runs L times slower than its maximal rate. Besides 
the resource waste, this renders multi-coset sampling impractical in wideband applications where /nyq 
is higher (typically by orders of magnitude) than the rate r of available devices. 

One recent paper [4] developed a nonconventional ADC design for wideband applications by means 
of high-rate optical devices. The hybrid optic-electronic system allows sampling at rate l/(LrNYQ) with 
minimal attenuation to higher frequencies (up to /nyq/2). Unfortunately, at present, this performance 
caimot be achieved with purely electronic technology. Thus, for wideband appUcations that caimot afford 
the size or expense of an optical system, multi-coset sampling becomes impractical. 

Another limitation of multi-coset sampling, which also exists in the optical implementation, is main- 
taining accurate time delays between the ADCs of different cosets. Any uncertainty in these delays 
hobbles the recovery from the sampled sequences. 

Before describing the way our proposed sampling stage overcomes these limitations, we briefly review 
the mechanism underlying the bhnd reconstruction of [3]. 

C. IMV System 

Let x(A) = {x(A) : A G A} be a collection of n-dimensional vectors indexed by a fixed set A 
that may be infinite. The support of a vector is the set supp(v) = {i | Vj / 0}, and we define 
supp(x(A)) = U;^ supp(x(A)). We will assume that the vectors in x(A) are jointly i^-sparse in the 
sense that | supp(x(A))| < In words, the nonzero entries of each vector x(A) he within a set of at 
most K indices. 

Let A be an m X n matrix with m <n, and consider a parameterized family of Unear systems 

y(A) = Ax(A), A G A. (2) 

When the support S = supp(x(A)) is known, recovering x(A) from the known vector set y(A) = {y(A) : 
A G A} is possible if the submatrix As, consisting of the columns of A indicated by S, has full column 
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Reconstruct joint support 



y(A) 



Construct a frame 
V for y(A) 



Solve V = AU for 
sparsest matrix U 



5 = Usupp(Ui) 



Fig. 2: Recovery of the joint support S = supp(x(A)). 



rank. In this case, 

X5(A) = (A5)V(A) (3a) 

Xi(A) = 0, i^S (3b) 

where X5(A) contains only the entries of x indexed by S, where A;^ denotes the conjugate transpose of 
As and where (A5)1^ = (A^A5)^^A^ is the Moore-Penrose pseudoinverse. For unknown support S, 
(2) is still invertible if K = jS*] is known, and every set of 2K columns from A is linearly independent 
[6]-[8]. In general, solving (2) for x(A) is NP-hard because it may require a combinatorial search. 
Nevertheless, recent advances in compressive sampling and sparse approximation delineate situations 
where polynomial-time recovery algorithms correctly identify supp(x(A)) for finite A. This challenge is 
sometimes referred to as a multiple measurement vectors (MMV) problem [8]-[13]. 

Recovering a multiband signal x{t) from a set of multi-coset samples can be reduced to a certain 
infinite measurement vectors (IMV) problem (where A is infinite). When the band locations are known, 
the support set S is determined and reconstruction can be performed via (3) [1], [2]. In a blind scenario, 
the support of the unknown vectors x(A) can be recovered in two steps [3], [6]. First, construct a (finite) 
frame V for y(A). Then, find the (unique) solution U to the MMV system V = AU that has the fewest 
nonzero rows. It holds that supp(x(A)) is the set S = Ui supp(Ui), where the union occurs over columns 
of U. Fig. 2 summarizes these recovery steps. 

In the next section, we describe and analyze the proposed sampling system. In contrast to the multi- 
coset strategy, our system uses standard low-rate ADCs. We match the analog input of the ADCs to 
their maximal rate. The system also avoids time offsets between devices. As in multi-coset sampUng, the 
sampling sequences generated by our system are related to x{t) via an IMV system, different from the 
one which is based on the sequences (1). Consequently, the recovery of x{t) can be performed via the 
steps described in Fig. 2 and (3). 
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Fig. 3: Description of a practical sampling stage for multiband signals. 



III. Efficient Sampling 

A. Description 

Let us present the proposed system in more detail. A block diagram appears in Fig. 3. We discuss the 
choice of system parameters in the sequel. 

The signal x{t) enters m channels simultaneously. In the zth channel, x{t) is multiphed by a mixing 
function Pi{t), which is a Tp-periodic, piecewise constant function that alternates between the levels ±1 
for each of M equal time intervals. Formally, 

Pi{t) = aik, k^<t<{k + l)^, 0<k<M-l, (4) 

with G {+1, —1}, and Pi{t + nTp) = pi{t) for every n G Z. 

After mixing, the output is converted to digital using the standard approach. In each channel, the signal 
spectrum is truncated by a lowpass filter with cutoff l/{2Ts) and the filtered signal is sampled at rate 

l/Ts. 

Note that the cutoff and the sampling rate match, and each channel operates independently. Since 
there are m channels, the average sampling rate is m/Tg samples/sec. A further advantage of this type 
of system is that samples are produced at a constant rate, so they may be fed to a digital processor 
operating at the same frequency, whereas multi-coset samphng requires an additional hardware buffer to 
synchronize the nonuniform sequences. 
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B. Analysis 

To ease exposition we choose an odd M, T = M//nyq, and Tg = Tp = T. These choices are relaxed 
in the sequel. Consider the i\h channel. Since pi{t) is periodic, it has a Fourier expansion 

oo 

Pi{t)= Q,e^'^"*, (5) 



where 



^ £ Pi(t)e-^¥n^dt (6) 

1 i-T/M 

= ;F / E (7) 

= 77; E "ifee-^t"'^ / e-^'^"*di. (8) 



Evaluating the integral gives 



1 . A 1 _ e-J-'o- 
Cm = E ^ , (9) 

where = 2it/M and Ci„ = Ci,_n. Expressing the Fourier transform Pi{f) in terms of the Fourier 
series coefficients Cj„ leads to 



/OO 



(10) 



with S{t) denoting the Dirac delta function. The analog multiplication Xi{t) = x{t)pi{t) translates to 
convolution in the frequency domain, 

00 

Xiif) = X{f) * Pi{f) = E Cin^ (/ - ^) • (1 1) 

n=— 00 

Therefore, Xi{f) is a linear combination of shifted copies of X{f). 

Filtering Xi{f) by H{f), whose frequency response is an ideal rect function in the interval JFq = 
[-l/(2r),l/(2r)], results in 

no 

Aif) = H{f)Xi{f) = E c^„x(/-^), feTo, (12) 

n=—no 

where no is the smallest integer satisfying 

2no + 1 > T/nyq. (13) 



8 



Under the choices above, no = {M — l)/2. The discrete-time Fourier transform of ai[n] is 



oo 



(14) 



n=—oo 



no 



0- 



(15) 



n=— no 



Substituting (8) in (15) leads to the system 



y(/) = (SF)(Dx(/)), 



(16) 



where yi(/) = Ai{e^'^'^^'^), l<i<m. Sis mmxM matrix whose ikth entry Sik = aik. The M x M 
matrix F is a certain cycUc columns shift of the discrete Fourier transform matrix of order M. The 
M-square diagonal D scales 



according to the last term in (8). Since D has non-zero diagonal entries, it can be absorbed into x(/) 
while keeping supp(x(jrQ)) = supp(Dx(jro)). Thus, (16) is an IMV system with SF replacing A of 



We now explain the choices assumed in the begiiming of the section. 

1) Tg = Tp justifies the imphcation (14)-(15). For Tg < Tp, (15) does not contain the entire information 
about X{f). Thus, this choice should be avoided. In contrast, the selection Tg > Tp has the following 
benefit for hardware implementation. Suppose Tg = kTp for some integer k > 1. Then, it can be 
verified that under technical conditions on Oik, each sequence ai[n] corresponds to k channels which 
are designed with Tg = Tp. Consequently, the number of mixers, filters and ADCs is reduced by a 
factor of k. This choice requires an ADC with r = kTp and is thus limited by the sampling rates 
of available devices. 

2) Using T = M//nyq ensures that the right hand side of (13) is an integer. Other choices are 
possible, but imply a sampUng rate that is higher than the minimum. 

3) Odd M simpUfy the exposition. However, a unique solution for the IMV (16) can also be guaranteed 
for even M. An even M also requires slight modifications for the reconstruction algorithms of [3] 
as will be detailed in the journal version of this paper. 

C. Parameter Selection and Stable Recovery 

The following theorem suggests a parameter selection for which the infinite sequences ai[n\,l < i < m 
match a unique x{t) G M.. When the band locations are known, the same selection works with half as 



Xi(/) = X(/ + (z-no-l)/r) 



(17) 



(2). 
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x(/) = 



HA 

Odd M 



x(/) = 



n 



HA 

Even M 

Tto — "2" 



Fig. 4: The relation between the Fourier transform X{f) and the unknown vector set x(/), eq. (17). 



many sampling channels. Thus, the system appearing in Fig. 3 can also replace the multi-coset stage of 
[2]- 

Theorem I (Uniqueness): Let a;(t) G be a multiband signal and assume the choices T = M//nyq 
for an integer M (not necessarily odd) and Tp = Tg = T. If 

1) M<fnYQ/B, 

2) m> 2N for non-blind reconstruction or m > 4N for blind, 

3) S = {aik} is such that every AN columns are linearly independent, 

then, for every / G To, the vector x(/) is the unique 2iV-sparse solution of (16). In addition, under these 
choices x(jFo) is jointly 4A^-sparse. 

Proof. The proof goes along the line of [3]. The relation (17) can be thought of slicing the spectrum 
X{f) into pieces of length 1/T and then rearranging them in a vector form x(/). Fig. 4 visuahzes this 
relation for even and odd M. 

The choice M < /mq/B ensures that 1/T > B and thus every band can contribute only a single 
non-zero value to x(/). As a consequence, x(/) is 2Ar-sparse for every / G To- In addition, this choice 
of M and the continuity of the bands guarantee that each band can occupy two spectrum pieces at the 
most. Therefore, when aggregating the frequencies to compute S = supp(x(jro)), we have \S\ < AN. 

In the non-bUnd setting, the band locations imply the support set S. The other two conditions on m, S 
ensure the existence of {As)^, and thus (3) provides the uniqueness of x(/). 
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In blind recovery, S is unknown, and the following CS result is used to ensure the uniqueness. A 
K-sparse vector x is the unique solution of y = Ax if every 2K columns of A are linearly independent 
[7]. Clearly, this condition translates to m > 4N and the condition on S of the theorem. ■ 

The parameter selection of Theorem 1 guarantees an average sampUng rate m/T > 4NB. Depending 
on whether /nyq/^ is an integer, this selection allows to achieve the minimal rate when taking the 
extreme values for m, M. Note that each x(/) is 2A?^-sparse, while x(jro) is jointly 4Ar-sparse under the 
parameter selection of the theorem. As detailed in [3], this factor requires doubUng m in order to use 
Fig. 2 and (3). Gaining back this factor at the expense of a higher recovery complexity is also described 
in [3]. 

Verifying that a set of signs {aik] satisfies the requirement of the theorem is computationally difficult 
because one must check the rank of every set of AN columns from S. It is known that a random choice 
of signs will work, except with probabihty exponentially small in M [14]. 

In fact, recent work on compressive sampUng shows that a random choice of signs ensures that signal 
acquisition is stable [9]. A matrix A is said to have the restricted isometry property (RIP) of order K, 
if there exists < (J/^- < 1 such that 

< IIAxf <(l+5;^)||xf (18) 

for every i^-sparse vector x [9]. When A = SF satisfies the RIP of order AN, then the matrices A5 
and {A.s)^ are well conditioned for every possible frequency subset C JTq with l^l < 2N. It was 
proved in [15] that basis pursuit can recover the X-sparse solution x of y = Ax for an underdetermined 
A, if A has < ^/^ - 1- The mean squared error of the recovery in the presence of noise or model 
mismatch was also shown to be bounded under the same condition. Similar conditions were shown to hold 
for other recovery algorithms. In particular, [16] proved a similar argument for a mixed ^2/^1 program, 
for the MMV setting. Thus, if A = SF has the RIP of order AN, then it implies the stability of the 
recovery using Fig. 2, when the mixed-norm program is utilized to solve the sparsest solution of the 
MMV V = AU. 

It remains to quantify when stable recovery is possible for specific choices of m and M. Let F be 
an M X M unitary matrix (such as F in (16)), and suppose that S is an m x M random matrix whose 
entries are equally likely to be ±.\/y/m. The RIP of order K holds with high probability for the matrix 

A = SF when m > CKlog{M/K), where C is a positive constant independent of everything [17]. The 
log factor is necessary [18]. In practice, we empirically evaluate the stability of the system since the RIP 
cannot be verified computationally. 
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IV. Numerical Evaluation 



To evaluate the empirical peri'ormance of the proposed system (see Fig. 3), we can simulate the action 
of the system on test signals contaminated with white Gaussian noise. To recover the signals from the 
sequences of samples, we apply the reduction from an IMV system to an MMV system, as described in 
Fig. 2. We solve the resulting MMV systems using simultaneous orthogonal matching pursuit [11], [12]. 

More precisely, we evaluate the performance on 100 noisy test signals of the form x{t) + w{t), where 
cc is a multiband signal and w is a white Gaussian noise process. The multiband signals consist of iV = 3 
pairs of bands, each of width B = 40 MHz, constructed using the formula 



where sinc(x) = sin(7rx) / (ttx). The energy coefficients are fixed Ei = {1, 2, 3}, whereas for every signal 
the carriers fi are chosen uniformly at random in [— /nyq/2, /nyq/2] with /nyq = 10 GHz. To represent 
the continuous signals in simulation, we place a dense grid of 4000 equispaced points in the time interval 
[— 200//nyq, 200//^fYQ]. The Gaussian noise is added and scaled so that the test signal has the desired 
signal-to-noise ratio (SNR), where we define the SNR to be 101og(||x||/||'u;||). 

We simulate the proposed system with m = 51 channels, where each mixing function pi{t) alternated 
sign at most M = 51 times. The sampling rate parameters are chosen so that Tg = Tp = M//nyq- 
Each sign Oik is chosen uniformly at random and fixed for the duration of the experiment. To simulate 
the analog lowpass filter, we use a 50-tap digital FIR filter, designed with the MATLAB command 
h=firl(50,l/M). The output of the filter is decimated to produce the sampled sequences ai[n]. 

The input signal is reconstructed from rh < m channels. We follow the procedure described in Fig. 2 
to obtain an estimated support set S. When S = S, the true support set, we declare that the system has 
recovered the signal. Fig. 5 reports the percentage of recoveries for various numbers fh of channels and 
various SNRs. 

To construct the frame V, we begin by computing the mP' values Qik = ^„ ai[n]afc[?7.]. We then 
perform the eigenvalue decomposition Q = VV^ and then discard the eigenvectors of the noise space 



N 




i=l 



[3]. 



V. Conclusions 

We developed an efficient sampling stage for analog multiband signals. In the proposed system, analog 
mixers and standard ADCs replace impractical nonuniform sampUng of multi-coset strategy. Analog 
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Fig. 5: Image intensity represents percentage of correct support set recovery S = S, for reconstruction from different 
number of sampling sequences rh and under several SNR levels. 

mixers for wideband applications is an existing RF technology, though selecting the exact devices may 
require an expertise in analog design. 

The proposed system has a set of parameters, which determines the signal, if selected according to the 
conditions we derived. Analyzing our system in the frequency domain lead to an IMV system, which 
allows to use existing reconstruction stages with only minor modifications. In addition, based on the IMV 
system and recent works in the CS literature, we deduce the rate requirements for stable bUnd recovery, 
which in general is higher than the rate required to determine the signal from its samples. 

A preliminary computer evaluation of our system shows a promise for stable blind recovery from 
sub-Nyquist sampling rate, although further work is required to quantify the optimal working point in 
the trade-off between sampling rate, blindness, and practical implementation. 
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